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ABSTRACT 

Water and simple organic molecular ices dominate the mass of solid materials available for plan- 
etesimal and planet formation beyond the water snow line. Here we analyze ALMA long baseline 
2.9, 1.3 and 0.87mm continuum images of the young star HL Tan, and suggest that the emission 
dips observed are due to rapid pebble growth around the condensation fronts of abundant volatile 
species. Specifically, we show that the prominent innermost dip at 13 AU is spatially resolved in the 
0.87mm image, and its center radius is coincident with the expected mid-plane condensation front of 
water ice. In addition, two other prominent dips, at distances of 32 and 63 AU, cover the mid-plane 
condensation fronts of pure ammonia or ammonia hydrates and clathrate hydrates (especially with 
CO and N 2 ) formed from amorphous water ice. The spectral index map of HL Tan between 1.3 and 
0.87mm shows that the flux ratios inside the dips are statistically larger than those of nearby regions 
in the disk. This variation can be explained by a model with two dust populations, where most of 
solid mass resides in a component that has grown into decimeter size scales inside the dips. Such 
growth is in accord with recent numerical simulations of volatile condensation, dust coagulation and 
settling. 

Subject headings: astrochemistry - stars: premain-sequence - stars: individual (HL Tau) - proto¬ 
planetary disks - planets and satellites: composition 


1. INTRODUCTION 

The frost lines of abundant ices have long been thought 
to be important for planet formation, in part because 
the extra material provided by condensation can signifi¬ 
cantly e nhance the local mass surface density of solids in 
the d isk (|Hayashill 19811 : [Stevenson fc Luninelil988HDeschl 
I2007f) . This in turn significantly shortens the time scale 
needed to form pr otoplanets in core-accretion models 
([Pollack et al.l[r99^ . Furthermore, the bulk composition 
of terrestrial planets and the cores of gas or ice giants is 
determined by that of the planetesimals from which they 
are built. Because volatility is highly species specific, 
it is the condensation fronts of the principal solid-vapor 
reservoirs in the mid-plane of protoplanetary disks that 
set the rough boundaries of planetesimals with different 
compositions (jOberg et al.ll201lll . 

Recent experiments and numerical simulations of dust 
coagulation and settling suggest that condensation fronts 
may play a crucial role in planetesimal formation via 
two aspects: First, the growth efficiency of dust grains is 
composition dependent. Lab experiments show that icy 
aggregates are significantly more ‘sticky’ and resistant 
to compaction than silicate aggregates (iGiittler et al.l 
120101: iSeizinger fc Klevll20I3l: iKelling et al.l 1201^ The 
second, perhaps more important, aspect is that the 
enhanced local surface mass density near snow lines 
may produce pressure bumps, high dust-to-gas mass 
ratios, or viscosity gradients, thus triggering insta- 
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bilities that form planetesimals (iKretke fc LinI l2007t 
Brauer et al.l l200§: iRos &: JohansenI l2013t Bitsch et al l 


2014 Drazkowska fc Dullemondll20l4iFlock et al.ll20I5ll . 

Despite the wealth of theoretical predictions concern¬ 
ing condensation fronts in disks, no direct observations of 
fast dust growth near such fronts have been carried out. 
The small angular separation and emitting areas of such 
zones, even in the nearest protoplanetary disks, make 
such observations challenging. ALMA, with its unprece¬ 
dented spatial resolution and sensitivity, is changing this 
situation dramatically. Young, actively accreting objects 
are good initial targets because their luminosities and 
disk temperature profiles move the condensation fronts 
out to larger, resolvable, d istances ([Kennedy fc Kenvoni 
l2nnaiBlake fc Berginl[201^ . 


HL Tau, a young star in the Taurus molecular cloud 
(~10®yr), was recently observed by ALMA at 2.9, 1.3 
and 0.87 mm with a spa tial resolution as good as 3 AU 
([Partnership et ahi I2015L APR henceforth). Here we 
propose that the dark concentric rings observed in the 
(sub)mm interferometry of HL Tau are due to fast peb¬ 
ble growth near condensation fronts. We first character¬ 
ize the location and width of surface brightness dips in 
the ALMA images before comparing their location with 
the expected condensation fronts in the mid-plane of the 
disk. We then examine the dust properties in the dips to 
determine if they are consistent with pebble growth. 


2. OBSERVATIONS 

The ALMA Band 3, 6 and 7 (or wavelengths of 2.9, 
1.3 and 0.87mm, respectively) HL Tau observations were 
carried out in 2014, as part of its long-baseline com¬ 
missioning and science verification (SV) program. With 
longest baselines of 15.2 km, the HL Tau observations 
achieved exceptional spatial resolution, with beam sizes 
of 11.9x8.6AU, 4.9x3.1AU and 4.2x2.7AU at 2.9, 1.3 
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and 0.87mm (for an adopted distance of 140pc). More 
detailed descriptions of the data and calibration can be 
found in APR, who demonstrated the HL Tau contin¬ 
uum images are ~axisymmetric and well fit with a series 
of elliptic rings. Therefore, a first order approximation 
of the HL Tau surface brightness is a single parameter 
function of the distance from the star. 

For our purposes, it is much easier to work on depro- 
jected, circular images. We thus start wi th the calibrated 
meas urements, and use CASA v4.3 (|McMullin et al.l 
1200711 to deproject the visibility data and generate syn¬ 
thesized images for all three bands. Assuming the peak 
continuum emission represents the location of the central 
star, we first shift the phase center of each band to that 
of HL Tau (by Aq, = —14 mas, = 196 mas). We then 
deproject the visibility data using the best-fit inclination 
angle i = 46.72° and PA = 138.02° from APR. 

From the deprojected visibility data in bands 6 and 
7 we have also generated a spectral index (a, where 
I,y oc map. This result is computed using the CLEAN 
task with nterms=2, based on th e multi-frequencv decon¬ 
volution algorithm developed bv iRau fc CornweTj (j2011l l. 

3. CHARACTERING THE EMISSION DIPS 

We measure the surface brightness versus radial dis¬ 
tance in steps of 0.5 AU, and present the azimuthally 
averaged radial brightness distributions in Figure [TJ The 
most consistent features in all the bands are three emis¬ 
sion depressions, or dips, around ~13, 32 and 63 AU. We 
note that there are many subtle bumps in the surface 
brightness distribution, and we refer interested readers 
to a more detailed description by APR. Here we focus on 
the three most prominent, ^radially-symmetric dips. 

In order to constrain the dip widths and depths, 
a smooth fit to the surface brightness distribution is 
needed. We use a simple vertically isothermal model, 

/4r)=cos^xH,(T,)(l-e-^W"■^ (1) 

r, = To(r/[lAU])« (2) 

E(r) = Ec(r/rc)“^exp[-(r/rc)^“'^] (3) 

where is the surface brightness (Jy arcsec“^), i the 
inclination angle, k^, the continuum opacity in cm^ 
and E(r) the dust-l-gas mas s surface density distribution 
for a steady accretion disk (jPringlelll981l : [Andrews et al.l 
1200911 . The cosi factor in eq. (1) accounts for the effects 
of the deprojection on the synthesized beam. 

Since our goal is to produce a smooth fit to the surface 
brightness distribution rather than a realistic physical 
model, we us e k„ = 0.01 cm^ g~^(t//230fG Hz]) through¬ 
out the disk (jOssenkonf fc HennindllOQ^ l. 

The results of the best-fit models are plotted as dashed 
lines in Figure [T] We then divide the observed Iv{r) with 
the best-fit model in each band to constrain the shape of 
the dips, see Figure [1] (b). The locations are consistent 
with those of the dark rings reported in APR. The 13 AU 
and 32 AU intensity dips appear to be highly symmetric. 
Gaussian fits to the normalized radial brightness profile 
are listed in Table[TJ The 63 AU dip is asymmetric, so we 
only list the radius of the minimum and the dip FWHM. 
The FWHM of the 13 and 32 AU dips are 5 —7 AU, about 
2x the minor axis of the synthesized beams at 1.3 and 
0.87mm, indicating that the dips are spatially resolved. 
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Figure 1. (a) Normalized radial surface brightness distributions 
of HL Tau at 0.87 (blue), 1.3 (red), and 2.9mm (black). The 
solid lines depict the observations, dashed lines the best-fit models 
(described in §3). (b) The observation-to-model surface brightness 
ratio. The three dotted vertical lines indicate the dip minima. 


Table 1 

Parameters of the dips, normalized to the local continuum 


band 

ind 

Center [AU] 

FWHM [AU] 

Amplitude 

band 7 

1 

12.6±0.1 

7.2±0.3 

0.56±0.02 

3.0X5.5AU 

2 

32.2±0.4 

5.3±1.0 

0.26±0.05 


3 

63.5 

17.7 

0.44 

band 6 

1 

12.7±0.2 

7.4±0.3 

0.66±0.03 

3.3X6.6AU 

2 

32.2±0.3 

6.9±0.6 

0.41±0.04 


3 

64.0 

18.0 

0.52 

band 3 

1 

13.6±0.3 

10.4±0.4 

0.39±0.03 

9.1x15.1 AU 

2 

32.4±0.6 

11.1±1.5 

0.22±0.03 


3 

66.0 

16.8 

0.45 


The 63 AU dip is resolved at all bands, with a FWHM 
of 17-18 AU. The minimum surface intensity of the dips 
range from 37% to 78% that of the smooth continuum 
from the best-fit models. 

4. GONDENSATION FRONTS OF MAJOR 
VOLATILES IN PROTOPLANETARY DISKS 

To check if the dark rings/dips in the ALMA HL Tau 
images are correlated with condensation fronts, we calcu¬ 
late the condensation temperatures of major ice species 
under realistic pressure ranges and compare them with 
the expected disk mid-plane temperature distribution. 

Here we consider the major condensible carriers of 
these elements based on the composition of comets, 
which are believed to be the best representatives 
of primordial icy mater ials in the Solar Nebula 
(|Mumma fc Charnlevll201ll) . We include N 2 , one of the 
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Table 2 

Condensation temperatures of the major volatiles in disks 


Species 


Eb 

Cometary Abundance 

Ref 


(Sy 

(K) 

% of H2O 


H20 

128-155 

5165 

100 

1,5 

CO 

23-28 

890 

0.4-30 

1,5 

C02 

60-72 

2605 

2-30 

1,5 

CH4 

26-32 

1000 

0.4-1.6 

2,5 

CH30H 

94-110 

4355 

0.2-7 

1,5 

N2 

12-15 

520 


2,5 

NH3 

74-86 

2965 

0.2-1.4 

1,5 

HCN 

100-120 

4170 

0.1-0.6 

3,5 

H2S 

45-52 

1800 

0.1-0.6 

4,5 

NH3-H20 

78-81 



6 

H2S* 

77-80 



6 

CHJ 

55-56 (69-72) 



6,7 

CO* 

45-46 (58-61) 



6,7 

N2* 

41-43 (55-57) 
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Note. — a. Condensation temperature ranges for ices correspond¬ 
ing to gas number densities of 10^ - lO^^cm”^, suitable for disk mid¬ 
planes. * Condensation temperatures of clathrate hydrate formed from 
hexagonal ice or hydrate under gas number densities of 10^^- lO^^cm”^. 
The values in pare ntheses are for clathrates f or med f rom am orphous ice. 
References: (1) IMartin-Domenech et al.l ll2014lh (2) ILuna et al.l 

ll20^^, (3) ISandford &: Allamandola| l| 1993l) . ( 4)JHa£eg ^wafcIierbs3 
(5) |M umma^^^h arnlevr ~| |201in , ^6) flr^T^dl |200^h (7) 
[Lumn^ fc StevensonTT^SSir 

major carriers of N as predicte d by chemical models of 
disks (iSchwarz fc Berginl l2014f) . As a homonuclear di¬ 
atomic, N 2 cannot be directly measured in cometary co¬ 
mae. 

For “pure” (single component) ices we define the con¬ 
densation temperature of species i as that where the ther¬ 
mal desorption and accretion rates are equal, i.e., 

^ice ^ ^desorp ^gas ^ ^accr 7 (^) 

where ^deso rp ^accr Calculated using the treatment 
outlined in iWoitke et al.l (I2009D and iWalsh et al.l (|2010D . 

Besides pure ices, we also consider clathrate hydrates, 
special forms of crystalline water ice in which gaseous 
molecules can be trapped inside lattice cages of the ice. 
At typical mid-plane disk pressures, abundant volatiles 
such as methane, H 2 S, CO and N 2 are expected to be 
trapped in the form of clathrate hydrates before they 
condense as pure ice species, given a sufficiently large 
water ice surface area to which the vapor has access 
(jLunine fc Stevensonlll985l l. This interesting feature of 
clathrate hydrates has been employed to explain the low 
N/0 elemental ratio observed i n comets as com pared to 
that of the Sun’s photosphere (liro et al.ll2003D. and th e 
existence of methane in Titan’s atmosphere (iLe'^ ilHnD- 

A summary of the condensation temperatures of pure 
ice species from such calculations is shown in Table [H 
along with th e condensation temp e ratur es of clathrate 
hydra tes from iLunine fc StevensonI (I1985D and liro et al.l 
(1200311 . Our conde nsation tempe r ature s are consistent 
with the results of iPollack et al.l (1199111 over the same 
pressure range. Table [5] shows that most of the clathrate 
hydrates can be created at higher temperatures than 
their pure ice condensates. The exceptions are NH 3 , 
which has a very similar temperature for its hydrate and 
pure ice forms, and CO 2 , which condenses as a pure frost 
at higher temperatures than for which the clathrate is 
stable. 



of HL Tau. The shaded areas show the clathrate-hydrate (abbrevi¬ 
ated as CL in the legend) condensation temperature ranges given 
in Table [2] the vertical lines the mean condensation front radii. 
The thick solid curve is the brightness temperature distribution of 
the 0.87 mm image, the thin curve the mid-plane temperature from 
MFH. A black diamond denotes the condensation temperature of 
pure CO2 frost, which is close to edge of the 63 AU dip. 

In Figure[2]we compare the dip radii with the expected 
disk condensation front locations. We also plot the ob¬ 
served brightness profile, which provides a rigorous lower 
bound to the mid-plane dust temperature. This pro¬ 
file would yield a water ice line at <4-5 AU, unresolved 
by the ALMA SV data, with the NH 3 /ammonia-hydrate 
condensation front near the 12 AU dip. 

The actual physical temperature in this embedded disk 
should be larger, however, especially in the case of opti¬ 
cally thin dust emission or significant grain coagulation. 
As a more detailed estimate of the HL Tau dust tem - 
perature structure we follow IMen’shchikov et al.l ()1999[ ) 
(MFH henceforth), who used two-dimensional radiative 
transfer models to quantitatively match the available 
spectral energy distribution, intensity, and linear po¬ 
larization data on HL Tau from near IR to (sub)mm 
wavelengths. Their best-fit model results in Tmid = 665 
(r/AU)“° ®K for the mid-plane temperature distribu¬ 
tion. This is greater than that estimated from the ob¬ 
served brightness distribution in the ALMA image, at all 
radii. Interestingly, as Figure [2] shows the dips overlap 
nicely with the condensation fronts of the most abundant 
volatiles using the MFH profile: pure water condenses 
around the 13 AU dip, while pure NH 3 or ammonia (and 
hydrogen sulfide) hydrates condense around the 32 AU 
dip. At further distances, the condensation front of pure 
CO 2 and the onset of the CO and N 2 clathrate hydrate 
stability fields, from an amorphous water ice seed occur¬ 
ring near the 63 AU dip. 

5. DUST PROPERTIES INSIDE THE DIPS 

We have shown above that the most prominent dips in 
the HL Tau images are remarkably close to the expected 
condensation fronts of abundant volatiles. We now in¬ 
vestigate the dust properties across these radii. 

At (sub)mm wavelengths, the spectral index a of the 
surface brightness (where oc iz“) is a widely-used ob- 
servab le to characterize d ust properties in protoplanetary 
disks (|Natta et al.ll200^ . Eor optically thin dust emis- 
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Figure 3. Panels (a)-(d): crosscuts of the spectral index a, measured between bands 6 and 7, at four angles (0 starts from the west). The 
grey regions depict la uncertainty intervals. Panel (e): An azimuthal average of the continuum spectral index o. The grey region again 
shows the la uncertainty range. Panel (f): Spectral index o as a function of average dust size, a, in the large dust aggregate population. 


sion, a can be expressed as a=2+/3, where k^, oc {v/vq)^ ■ 
The value of /3 is sensitive to the maximum dust size, 
Umax, for a given a size distribution, n(a) oc a'^. If the 
emission is optically thick, however, a=2 (the blackbody 
limit) and information on the dust size is lost. 

APR provided a cross-cut of a along the major axis of 
the as-observed image; but it is not clear if this cross-cut 
is representative of the full disk. In Figure [31 we show 
the radial distribution of a, measured between 1.3 and 
0.87mm, in crosscuts at 0 = 0, 45, 90 and 135° (start¬ 
ing from west) along with the azimuthally averaged a{r). 
Despite large azimuthal variations, there is a clear gen¬ 
eral pattern: inside of 48 AU a^2, except for increases 
to values of ^2.6 and ^2.4 at the 13 AU and 32 AU dips, 
respectively. Outside 48 AU, a gradually increases from 
2 to 2.7 and then becomes relatively flat, except for two 
shallow decreases around 70 and 80 AU near local max¬ 
ima in the surface brightness. 

One way to explain a~2 invokes optically thick dust 
emission in the inner disk, while that in first two 
dips is optically thin with a=/3-|-2. This explana¬ 
tion faces two difficulties. First, a massive disk is 
needed to create optically thick mm-continuum with 
modest grain growth. Given a typical continuum opac¬ 
ity of Ki,=0.01 cm^ g“^(z//[230GIIz]), a surface density of 
Sgas=100 g cm“^ at 48 AU is needed to reach = 

1 - a value 20 x larger than that for the MMSN (jHavashil 
Il981fl . The second issue is the large discrepancy be¬ 
tween the temperature derived from SED fits (the MFH 
model) and the observed brightness temperature (Fig¬ 
ure HI). These should be nearly equal if the mm-dust 
emission is optically thick. 

An alternative explanation for a^2 that mitigates 
these difficulties is a scenario where the dust emission 


for R <48 AU is effectively optically thin, but with /3~ 0. 
This can occur when dust grains grow to sizes (a) larger 
than the observational wavelength (A). Under such con¬ 
ditions, the absorption cross section simply equals the 
particle geometric cross section - and is thus indepen¬ 
dent of wavelength. Such significant dust growth is likely 
in the dense inner disk, considering that grains beyond 
48 AU have grown to >mm sizes (a=2.77±0.13 (APR)). 
This dust growth scenario, however, also faces a problem: 
for a single dust size-distribution, a > 2 within the 13/32 
AU dips demands smaller dust (and reduced mass sur¬ 
face densities) at these radii. It is difficult to understand 
why the disk inside of and, particularly beyond, the 13 
and 32 AU dips should have experienced more extensive 
dust growth than the dense inner regions of the disk. 

We suggest the apparent contradiction in the dust 
growth scenario can be solved by a model with two dust 
populations. This model differs from the commonly in¬ 
voked power law size distribution but is commonly en- 
countered in numerical simu l ations of dust growth (e.g. 
Dull emond fc DominikI 120051 iDrazkowska fc Dullemondl 
2014lb Our first dust population has intermediate 
growth, while the second population has grown to an 
average size for which a ^ A. Thus: 


^dustl,z/ — ^l( )^^ 

t'O 


/^dust2,2/ — ■ 


7ra 


A1 

4p a 


^dust,i/— (1 /)^dustl,i/ T /^dust2,i^ 


(5) 

( 6 ) 
(7) 


where a and / are the average dust size and mass frac¬ 
tion of the second dust population, and p is the physical 
density of the dust. 

A plausible order of magnitude estimate can be found 
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if we assume Kdusti.iv = 1 cm^ g“^(i//[230GHz]), or dust 
grain growth to mm size. For the second population, 
we assume p=lg cm“^, that is, ice dominated grains, 
and «;dust 2 ,i' =0.75(a/[cm])“^ cm^ g“^. We further adopt 
/=0.9, i.e., 90% of the dust mass in the second dust 
population. For modest optical depths in the first pop¬ 
ulation,the overall value of a between 1.3 and 0.87mm 
changes with d (Figure[3]f). For d near 1 cm, the cm-sized 
dust dominates the opacity and leads to 2. When d 
grows to decimeter-size, the mm-sized grains (population 
1 ) dominate the mass specific dust opacity, even though 
they contain a small fraction of the dust-l-ice mass, and 
a becomes significantly greater than 2. 

Thus, with significant dust aggregation in the second 
population, the flux and a behavior in the dips can be 
explained. With modest mass surface densities,the frac¬ 
tional area of a given column of the disk covered by 
the dust becomes <unity, and the brightness tempera¬ 
ture drops below the physical temperature even though 
individual dust aggregates are optically thick. Further 
(sub)mm to cm high resolution continuum images can 
distinguish between these two dust/disk structure sce¬ 
narios through more accurate constraints on the a{r) 
distribution, as can direct measurements of the vertical 
and radial gas temperature distributions via molecular 
emission. 

6 . DISCUSSION 

The analysis presented above demonstrates that the 
location of the three most prominent dips in the HL Tan 
(sub)mm interferometric images are coincident with the 
expected condensation fronts of the main volatiles, using 
the mid-plane temperature distribution derived from ob¬ 
servations. Further, the spectral index variation inside 
the dips can be explained by a bimodal dust size dis¬ 
tribution model without the need to invoke significant 
surface mass density depletions. 

From our simple order of magnitude model, we infer 
that most of the dust mass needs to reside in a popu¬ 
lation that has gro wn to decimeter s i ze sca les inside the 
13 and 32 AU dips. IRos fc Johansenl (I20I3I1 showed that 
near the water condensation front dust growth from mil¬ 
limeter to at least decimeter-sized pebbles is possible on 
a time scale of only 1000 years. This rapid growth is 
consist ent with young age of HL Tau, between O. I and 
IMyr (iBeckwith et al.l[l990l : IRobitaille et al.ll2007ll . 

To date, the numerical simulations for dust growth 
around condensation fronts have considered only water. 
It is unclear if pebble growth is sufficiently rapid around 
other condensation fronts to explain the HL Tau results. 
The extra solid material delivered by clathrate hydrates 
cannot significantly enhance the solid mass surface den¬ 
sity as does H 2 O at the water snowline, since each host 
molecule needs ~ six water molecules to form the cage 
structure. The condensation fronts of pure CO 2 and 
CO are potentially more important for enhancing mass 
surface density since they account for ~30% of water 
abundance in comets. Nevertheless, the model used by 
IRos fc Johansenl (1201311 can in principle to be applied to 
other condensation fronts. Clathrate hydrate formation, 
per haps driven by th e transient warming of amorphous 
ice (iBlake et al.llI99ill via accretion bursts, can alter the 
ice rheology. Changes to the sticking efficiency, poros¬ 
ity and compaction of dust-|-ice aggregates may thus be 


central to triggering rapid pebble growth in the outer 
disk. 

Once a critical decimeter-sized pebble population is 
formed, streaming instabi lities can drive the cr eation of 
>km-sized planetesimals ([Johansen et al.l[20I4ll . If fast 
pebble growth does preferentially occur around the con¬ 
densation fronts of abundant volatiles, this would sug¬ 
gest that snow lines regulate the formation and chemical 
composition of I-IOO km planetesimals, and ultimately 
the formation and bulk composition of planets. 

As the disk evolves the location of condensation fronts 
will shift inward, perhaps countered by the effects of 
episodic accretion . How this might affe ct dust evolu¬ 
tion is uncertain (iHubbard fc Ebelll20I4ll . Nonetheless, 
the possibility that disks such as that encircling HL Tau 
might be seeded with pebbles suggests that planetesimal 
formation might occur during early evolutionary stages. 
This is consistent with the cosmochemical record in our 
Solar system which shows that large and differentiated 
bodies had already begun to form <I Myr after the con¬ 
densation of calcium-aluminum-rich in clusions (CAIs) , 
the oldest minerals in the Solar system (|Qin et al.ll2008l : 
iKleine et al.ll2009l : IKruiier et al.ll20lll . 


This paper makes use of the following ALMA data 
sets: ADS/JAO.ALMA#20II.0.000I5.SV. ALMA is a 
partnership of ESC (representing its member states), 
NSF (USA) and NINS (Japan), together with NRC 
(Canada), NSC and ASIAA (Taiwan), and KASI (Re¬ 
public of Korea), in cooperation with the Republic of 
Chile. The Joint ALMA Observatory is operated by 
ESO, AUI/NRAO and NAOJ. The authors gratefully ac¬ 
knowledge support provided by the NSE Astronomy & 
Astrophysics, NSE INSPIRE (AST-1344133), and NASA 
Origins of Solar Systems grant programs. 
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